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ABSTRACT 



Distance relations in a locally inhomogeneous universe are expected 
to behave like the Dyer-Roeder solution on small angular scales and the 
Friedmann-Robertson- Walker solution on large angular scales. Within a simple 
compact clump model the transition between these asymptotic behaviors is 
demonstrated and quantified. The redshift dependent transition scale is of 
order a few arcseconds; this implies it should have little influence on large 
angular scale cosmological tests such as the volume-redshift relation but possibly 
significant effects on arcsecond angular diameter measurements of radio galaxies 
and AGNs. For example, at z = 2 on arcsecond scales a clumpy flat universe 
mimics the angular diameter distance of a smooth Q = 0.27 model. 



Subject headings: gravitational lensing, distance scale, cosmology: theory 



1. Introduction 

Light propagation through our universe is generally treated in two separate regimes: 
either within the smooth Friedmann-Robertson- Walker model or within the gravitational 
lensing model of density inhomogeneities (see, e.g., |Schneider, Ehlers, fc Falco 1992p . 



In the actual universe both local inhomogeneities and global homogeneity exist, so each 
acts as an asymptotic description under the appropriate conditions. Along with image 
amplifications, distortions, and time delays, the source distance-redshift relation is also 
affected by inhomogeneities. While many researchers discuss the distance relations in the 
presence of inhomogeneities, e.g. Futamase & Sasaki (1989), Watanabe & Tomita (1990), 
and how ignorance of dumpiness can affect determination of cosmological parameters by 
distance tests, e.g. Linder (1988a), Hadrovic & Binney (1997), the question of the nature 
of the transition region between the two regimes is not addressed. It is of importance to 
understand where this transition occurs so as to use the two limiting cases only where they 
are valid approximations. 
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For almost all astrophysical applications of light propagation geometric optics holds, 
where the curvature or inhomogeneity scales are much greater than the wavelength of the 
electromagnetic radiation, so the radiation can be treated as a beam or bundle of light rays. 
The focusing of the bundle due to the spacetime geometry determines the angular diameter 
distance r by 



d 2 r/d\ 2 = -(K+ \a\ 2 )r } (1) 

where A is the affine parameter measuring the path length, 1Z = (l/2)R flu k' J 'k u is the Ricci 
contribution of the gravitational focusing from matter within the beam, and o is the Weyl 
contribution of gravitational shear from matter outside the beam. Here is the Ricci 
tensor and fc M the photon four momentum. 

To obtain a universal, i.e. isotropic, distance relation, equation (]!]) is generally solved 
under the "average path" assumption (pyer fc Roeder 1973[) for an infinitesimal light 



beam, e.g. a single ray. This posits that the matter density along the line of sight is 
given by the global average density times a uniform smoothness parameter a, measuring 
the degree of small scale inhomogeneity or dumpiness, and that the shear vanishes from 
global homogeneity. This approximation substitutes for detailed knowledge of the locally 
inhomogeneous metric, or physical conditions, along the light path, which we generally lack. 
Given that ignorance some such effective model must be adopted, with the main caution 
being insurance that a "typical" path is indeed characteristic of the average. 

As the light beam subtends larger and larger solid angles, at some point the 
inhomogeneities should be smoothed over and we can legitimately calculate distances 
within the Friedmann- Robert son- Walker cosmology. The question arises how to connect 
these two asymptotic behaviors with a reasonable, preferably simple analytic model. (One 
can of course use numerical ray shooting within a specific pattern of inhomogeneities to 
compute the distance relations but they will be relations, not one single universal relation 
as desired.) Of particular observational interest is how large is the transition angle, as 
many cosmological tests such as the number-redshift and magnitude-redshift relations are 
sensitive to the precise distance measure. 



2. Smoothing 

The standard model for density inhomogeneities when calculating cosmological distance 
relations is the Dyer-Roeder prescription: distribute randomly a fraction 1 — a of the 
total matter density in clumps, point masses M with constant comoving number density 
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n . Light rays passing too close to a clump will be appreciably gravitationally influenced 
(spatially deflected, cross sectionally sheared, intensity amplified) and so will be recognized 
as something special and not treated as being on a typical line of sight. That leaves those 
paths avoiding clumps, but their neighborhoods only possess a matter density of ap where 
p is the global average matter density. Thus, an "average" path does not feel the full 
Friedmann-Robertson- Walker (FRW) density and the calculated angular diameter distance 
r(a,z) will differ from the FRW relation r FRW = r(l,z). In this case 71 = 4n(l + z) 2 ap 
where z is the source redshift. 

[One can also generalize the model to redshift dependent dumpiness a(z), modeling 
the effects of evolving inhomogeneity, clump mass, or number density ( [Linder 1988b[ ). Note 



that this brings up an interesting point. If we lived in a local inhomogeneity, either a 
bubble of different density or different dumpiness, then even though the volume averaged 
universe would appear Friedmann, our observed distance-redshift relation would not be 
the same as in that volume averaged universe. It does not asymptotically approach the 
Friedmann result even at distances much greater than the bubble size. I.e. the observer 
is in a special position and volume rather than angle (line of sight) averaging may give 
misleading results; one cannot invoke an ergodic homogeneity because our position as 
observers selects a unique location. Such an isotropic inhomogeneity would lead to an 
apparent conflict between cosmological parameters derived from dynamical quantities and 
those involving the distance-redshift relation (in preparation).] 

Given the Dyer-Roeder ansatz, one can see that as a light beam subtends a larger solid 
angle, the probability increases that it will include a clump. This raises the matter density 
within the beam and when the bundle is broad enough it will feel the full FRW density. To 
investigate this transition consider a beam with half angle 9 at the observer, propagating 
from a redshift z. Let each clump own a conical volume V c given by 

ndV c = l, (2) 

which extends from the observer out to the source or survey depth z with half angle 8 C . 
Here n is the proper number density; note the subscript c on the proper volume element dV 
denotes "clump" not "comoving". 

The probability that the light beam will include the clump is given by the Poisson 
process 



p(6) = l-e-^\ (3) 
where the optical depth r = V(9, z)/V(9 c (z)). That is, on average a beam covering V c 
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has roughly unity probability for including the clump (really 1 — e 1 because of statistical 
fluctuations in the number density). 

Although as the beam expands, i.e. 9 increases, the matter within the beam actually 
jumps in a step function, as the next step in our ansatz we simulate the global averaging 
by choosing an effective smoothed density p e ff = dp + p(l — a)p(9), where a is the global 
smoothness parameter, i.e. that felt by an infinitesimal ray. So rather than having no 
clump within the beam then suddenly all of it, we smoothly add its contribution according 
to the probability for inclusion p(9). This is the main approximation we make in order to 
talk about a universal distance function. It seems physically reasonable and is well behaved 
mathematically. 

The effective smoothness parameter is now 



a(9, z) — a + (1 — a) 
The volumes are given by 



l-exp{-V(9,z)/V[9 e (z)]} 



(4) 



V(9,z) = J dr p J g dA = J dr p J^dujr 2 ^)} 

= 2ir [ dr p [ d$$r 2 {a($)] = n9 2 g(z,a), (5) 

Jz JO 

where r p is the proper distance along the line of sight, A is the transverse area on the sky, 
and cu the angular area on the sky. 

Substituting equation (||) into @) yields 



a(6) 



a + (1 — a) 



1 - exp{-(6/9 c ) 2 [g(z,a)/g(z,a c )]} 



a apx (9) = a + (l-a)[l-exp{-(^ c ) 2 }], 



(6a) 
(6b) 



plotted in Figure 1, where a c = a( 



■a) 



We adopt Q — 1, which would cause 



the greatest deviation of the clumpy vs. smooth FRW model. Note that the approximation 
in (|6b|) (solid curve) has removed the need for recursion in determining a{9). Because the 
ratio of g's is insensitive to a the approximation is excellent: the fully recursed a differs 
from a apx by at most 11% for the extreme case of a = 0, z = 3 (and is more typically < 1%) 
- and this occurs in a region of parameter space such that the distance r(a) is never more 
than 1% from r(a apx ). 
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The resulting a(9) has the pleasing properties of monotonicity, simplicity, and the 
proper asymptotic behaviors: as 9 increases the parameter goes from a(9 <ti 9 C ) — a, to 
a(9 3> 9 C ) = 1, i.e. we have effectively introduced an averaging procedure that provides the 
Dyer-Roeder clumpy universe result at small angles, the FRW result at large angles, and 
defines the transition. 

To find the transition angle 9 C combine equations (0) and (|5|) to get 

1 = 2vm f C d9 9 [ drjl + z) 3 r 2 [a(9)] = n9 2 c n H^ f(z, a c ), (7) 

JO Jz 

for constant comoving clump density. So 



9 C = [nn H 3 f(z,a c )}- 1/2 . (8) 

Taking no = (3Hq/8tt)Q c M _1 where Ho = 100 h km s _1 Mpc -1 is the Hubble constant and 
Q c = (1 — a)Q the clump matter density in units of the FRW critical density, 



9 C = 1.3" {M/lO^M^h^il - a)- 1 / 2 /- 1 / 2 . (9) 

At low z or high smoothness 1-a < 1 the clump angle formally becomes large but this 
is unphysical given our interpretation of 9 C as a transition angle to FRW behavior. A more 
physical transition angle G is defined as that angle where r(a(9)) first differs from rpnw by 
less than 5%. 

Figure 2 plots both G and 9 C vs. z. For low z the distances r are insensitive to a (it 
enters at third order in an expansion of r in z; see [Linder 1988b|) so even for 9 = all r(a) 



are close to tfrw- For high redshifts 9 C and G formally diverge as a — > 1 due respectively to 
the scarcity of clumps and to the difference in high z asymptotic behavior between smooth 
and even slightly clumpy distances. The growth is extremely slow, however, 9 c (z ^> 1) ~ 25" 
for a = 0.999 and Q(z = 10 3 ) ~ 7" for a = 0.99. Realistically, a evolves due to structure 
formation such that it approaches unity sufficiently closely (e.g. > 0.998 for z < 10 3 ) 
that G — > at high redshifts. Thus, in general distant source observations with fields of 
view larger than a few arcseconds can legitimately be treated within the FRW model for 
distances. 

As a definite example, consider the volume-redshift cosmological test which involves 
the angular diameter distance (note the luminosity and proper motion distances are simply 
related to the angular diameter distance by factors of 1 + z). From equation (|5|) the 
volume-redshift relation is 
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V(9) = V FRW (9)[g(z,a)/g(z,l)}, (10) 

where a frw — 1 • For z < 1 the ratio of g factors is close to unity, and for z > 1 we observe 
at 9 3> 9 C (larger than arcsecond fields) and so a(9) ~ 1 and again V ~ Vfrw- Thus 
dumpiness effects are negligible in this case. 

The differential volume test is more sensitive since this does not possess the integration 
in equation (j|) that includes those very low redshifts where r is nearly independent of a (as 
discussed above). Then 



AV(9)/AV FRW (9) = 29- 2 f d$$r 2 [ a m/r 2 (l). (11) 

Jo 

The ratio is bounded between [r(a(9)) /r(l)] 2 and [r(0)/r(l)] 2 and can approach 1.25 at 
z = 1 as 9 — > if a(9) were unreasonably pushed to zero. The deviation from unity is less 
than 5% for all redshifts and dumpiness factors, however, when 9 > 12". In particular, this 
assures us that the usual observations can be analyzed properly using the FRW volume 



element: for example the Loh-Spillar (1986) fields were T x 10' at z < 0.75. (See |Omote fc 



Yoshida 1990| for how the clumpy luminosity distance affects the analysis of flux weighted 



counts. 



Clumpy and transition regimes can be important, however, for cosmo logical tests 
involving arcsecond scales. Possible applications include observations of radio galaxy lobes 
( |Kapahi 1989| , |Guerra fc Daly 1996[ ) and milliarcsecond observations of active galactic 



nuclei ( jGurvits 1994] ) for use in the angular diameter distance-redshift cosmological test. 



At a redshift of two, for example, the difference between the clumpy and smooth angular 
diameter distances for Q — 1 is 33%, and a clumpy flat model gives the same distance as a 
smooth model with Q = 0.27. 

Figure 3 shows this trade off between dumpiness in a Q = 1 model and low density in 
a smooth model. Thus a flat but clumpy universe could be misinterpreted through small 
angular scale observations as a lower density one. In general one can match a distance in a 
universe with a (a f > a,Q' < Q) model. The numbers on the plotted curves give the 
angular scale of the observations in arcseconds and show the transition in dumpiness or 
alternatively density miscalculation from the infinitesimal case (Dyer-Roeder) to the large 
scale average (FRW). 
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3. Shear 

One point remains of which we must be cautious. The TZ term in equation ([!]) was 
treated by using an effective smoothed density but we neglected the shear term \a\ 2 . On 
large scales we expect the shear to average to zero due to homogeneity but this is not 
ensured at smaller scales: if the light passes far from the clump then the shear on the 
bundle should be small, but not if it passes near the clump. 

We can analyze this in terms of an effective dumpiness defined through equation ([[]) by 

a eff = a- (2/3)\a\ 2 H 2 n-\l + z)' 5 , (12) 

where Q is the ratio of the total density to the critical density. Then the right hand 
side of equation (M can be written simply as 47r(l + z) 2 a e fjpr and the previous results 
hold with the substitution of a e ff for a. Note that shear always decreases a e // - brings 
it further from the FRW value of unity - and so increases the transition angle 8 C (more 
properly decreases the ratio 9/9 c for which a e ff is a given value). Also note that a e ff can 
be negative, which causes no mathematical worries [it merely makes the usual parameter 
P= (25 - 24a) 1 / 2 > 5]. 

Alternately, we can solve the full beam equation ([!]). Because the shear squared source 
term of a point mass dies off as (distance) -4 , we can approximate its behavior as localized in 
affine parameter, say between Ao and Ao + AA. Dividing the light propagation from source 
at A s to observer at into three regimes (0, Ao), (Ao, Ao + AA), (Ao + AA, A s ), we match 
the values of r and its first derivative at the boundaries. For the most extreme case, a = 
and Q = 1, the solutions to the distance relations are: r\ = A; r 2 = C sin |a|A + D cos \a\ A; 
r 3 = AX + B. Taylor expanding under the assumption A A <C A , | cr | ~ 1 yields 



r 3 (A) w A-|a| 2 A 2 AA(A/A -l) 

« n(A) - (4/25)|a| 2 A^ - 6 (l - y 5/2 )[l - (y /yf% (13) 



where y = 1 + z and A = (2/5) (1 — y 
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Thus the shear only has a significant effect on the angular diameter distance relation 
if the second term is nonnegligible. The maximum shear on the beam is given in order of 
magnitude by \a\ ~ Mb~ 2 where b is the impact parameter from the mass M. For b equal to 
the Einstein radius of the mass, \a\ ~ Ho(r s /riri s ) where the r's are respectively measured 
from observer to source, observer to lens mass, and lens to source. The maximum value of 
the ratio of the second to first terms in equation (|1^) is then (if ^o) 2 AA/A 1, meaning 
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that when concerned with distances we can neglect shear for beams passing outside the 
Einstein ring of the mass. 

So we argue that we can neglect the influence of shear on distances overall because 
1) narrow beams (9 <C Qe = bE/ri) would typically miss the Einstein ring because their 
probability for intersection (0 E /6 2 C ) is small (e.g. 9 2 E /9 2 C m Vt c z 2 <C 1 for z <C 1); 2) those few 
appreciably sheared light beams would be recognized as atypical and not used for distance 
measures; and 3) broad beams would have any shear effects diluted due to those portions of 
the beam that lie far from the clump and those symmetric (isotropic) about it. Therefore 
we do not expect shear to alter significantly the distance-redshift relation for a typical light 
bundle, except possibly at high redshifts where the transition scale 6 C is below an arcsecond 
(and so 9 C < 9 E ) and hence below the region of observational interest. 

4. Conclusion 

Using a simple toy model of an effective density distribution due to inhomogeneities 
one can derive a universal angular diameter distance relation applicable over all angular 
scales. It has the desired asymptotic properties: agreeing for small solid angles with the 
clumpy universe Dyer-Roeder distance, recreating for large solid angle observations the 
Friedmann-Robertson- Walker relation, and interpolating smoothly between them. The 
transition angle depends on the type of inhomogeneities but is estimated to be unlikely to 
exceed 10" for any cosmological test and should be of order 1" for any quantities involving 
the entire light propagation path between the source and observer. Still, this can have 
significant effects on such relations as the angular diameter-redshift relation for radio 
galaxies and AGNs, for example a true clumpy flat universe mimicking a smooth open 
model, and must be looked for in such cosmological tests. 

I thank Neal Katz for useful suggestions. This work was supported by NASA grants 
NAG5-3525 and NAG5-4064. 



- 9- 



REFERENCES 

Dyer, C. C, & Roeder, R. C. 1973, ApJ, 180, L31 
Futamase, T., & Sasaki, M. 1989, PRD, 40, 2502 

Guerra, E. J., & Daly, R. A. 1996, in Cygnus A: Study of a Radio Galaxy, ed. C. Carilli & 
D. Harris (Cambridge: Cambridge U. Press), 252 

Gurvits, L. I. 1994, ApJ, 425, 442 

Hadrovic, F., & Binney, J. 1997, MNRAS, submitted, ( |astro-ph/9708lT0|) 

Kapahi, V. K. 1989, AJ, 97, 1 

Linder, E. V. 1988a, A&A, 206, 175 

Linder, E. V. 1988b, A&A, 206, 190 

Loh, E. D., & Spillar, E. J. 1986, ApJ, 307, LI 

Omote, M., & Yoshida, H. 1990, ApJ, 361, 27 

Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses, (Berlin: Springer- Verlag) 
Watanabe, K., & Tomita, K. 1990, ApJ, 355, 1 



This preprint was prepared with the AAS IATgX macros v4.0. 



-10- 



Fig. 1. — The effective smoothness a is plotted vs. beam size 9 in units of the clump angle 
9 C . The lower curves are for a completely clumpy universe, a = 0, while the upper have 
a = 0.5. Solid lines illustrate a apx from (§b|), dotted and dashed lines use the full a from 
(|5a|), evaluated at survey depths of z — 1 and z = 3 respectively. 

Fig. 2.— The transition angle G (solid curves) and clump angle 9 C (dotted curves) in 
arcseconds are plotted vs. redshift. The pairs of curves are labeled near their intersections 
by the dumpiness a. At low redshifts the distance relation is close to the FRW behavior 
for all beam sizes while at high z it makes a transition to the clumpy behavior for angles 
smaller than 0. 

Fig. 3. — The dumpiness a needed for a flat universe to have the same angular diameter 
distance to redshift zasa smooth Q < 1 universe is plotted. In the conventional Dyer-Roeder 
model, for example, at z = 2 one could match r(z) for a = 0, Q = 1 by a FRW model with 
Q = 0.27. In the transition model of this paper, though, a is effectively a function of beam 
size. The numbers superposed on the curves give this size in arcseconds (for infinitesimal 
beam dumpiness a — 0). Thus observations with a 2" beam at z = 2 in a clumpy flat 
universe can be misinterpreted as belonging to a Q = 0.54 FRW universe. 
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